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For the effective-range function k cot 5, a pole approximation that involves a small number 

of parameters is derived on the basis of the Bargmann representation of the S matrix. The 

parameters of this representation, which have a clear physical meaning, are related to the pa- 
CD ■ 

(Sj . rameters of the Bargmann S matrix by simple equations. By using a polynomial least-squares 

fit to the function k cot 5 at low energies, the triplet low-energy parameters of neutron-proton 
scattering are obtained for the latest experimental data of Arndt et al. on phase shifts. The 
results are a t = 5.4030 fm, r t = 1.7494 fm, and t>2 = 0.163 fm 3 . With allowance for the values 
found for the low-energy scattering parameters and for the pole parameter, the pole approx- 
imation of the function k cot 5 provides an excellent description of the triplet phase shift for 

^ ■ neutron-proton scattering over a wide energy range (Ti ab < 1000 MeV), substantially im- 

o : 

. proving the description at low energies as well. For the experimental phase shifts of Arndt 

O ; 

et al., the triplet shape parameters v n of the effective-range expansion are obtained by using 

^— > ■ 

h ' the pole approximation. The description of the phase shift by means of the effective-range 

<-> : 

expansion featuring values found for the low-energy scattering parameters proves to be fair- 

■ 

>- ! ly accurate over a broad energy region extending to energy values approximately equal to 

•i-H . 

the energy at which this phase shift changes sign, this being indicative of a high accuracy 
and a considerable value of the effective-range expansion in describing experimental data 
on nucleon-nucleon scattering. The properties of the deuteron that were calculated by using 
various approximations of the effective-range function comply well with their experimental 
values. 
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1. INTRODUCTION 



The S matrix [1, 2] is a fundamental quantity in scattering theory. At a fixed value of 
the orbital angular momentum £, the S matrix is a function of the wave number k related 
to the energy E of two colliding particles in the cm. frame by the equation E = h 2 k 2 /2m , 
where m is the reduced mass of the system and h is the Plank constant. In the following, 
we restrict our consideration to the case of zero orbital angular momentum (we omit the 
index t = for the sake of simplicity) . The fact that a unified description of both scattering 
and bound states in two-particle systems can be constructed on the basis of the respective 
S matrix is an important consequence of the analytic properties of this matrix. In general, 
the analytic properties of the S matrix in the k plane are rather complicated; however, 
a number of important results can be obtained in considering specific physical systems by 
using an S matrix whose analytic properties are simple. In the k plane, the S matrix satisfies 
the relation [1] 

S* (k*) ■ S (k) = 1 . (1) 

At real k values, this relation coincides with the unitarity condition. 

In potential-scattering theory, the S matrix S (k) is expressed in terms of the Jost function 
F (k) as [3, 4] 

At complex values of k, the Jost function F (k) satisfies the condition 

F* = F (jfe) . (3) 

The S'-matrix property (1) immediately follows from this condition with allowance for (2). 
The behavior of the Jost function at high energies is determined by the relation 

lim F (k) — 1 . (4) 

The behavior of the S matrix at high energies directly follows from (2) and (4): 

lim S(k) = 1. (5) 

fe^oo 

In the case of elastic scattering, the S matrix can be expressed in terms of the phase shift 
S (k) as 

S (k) = e 2l<5(fe) . (6) 



Using relation (6) and considering that the S matrix possesses the symmetry property [1] 

S(-k) = S- 1 (k) , (7) 

one can easily see that, at real values of k, the phase shift 5 (k) is an odd function of k; that 
is, 

5 (—k) = -5 (k) . (8) 

From relations (5) and (6), one readily establishes the asymptotic behavior of the phase shift 
at high energies: 

lim 5 (k) = . (9) 

fe^oo 

In practice, it is very useful to introduce the effective-range function 

k cot 6 = ik S (*l + \ . (10) 
o [k) — 1 

In the physical region, the function k cot 5 is real by virtue of the unitarity of the S matrix, 
while, in the k plane, this function is analytic everywhere, with the exception of the points 
where S (k) = 1, i.e. the points where it has poles. The introduction of the function A; cot 5 
simplifies the investigation of the analytic properties of the S matrix. It follows from (8) 
that the function A; cot 5 is an even function of k. 

In studying nucleon-nucleon (NN) scattering at low energies, it is useful to employ the 
well-known effective-range expansion 

k cot 5 = -- + \r k 2 + v 2 k^ + v 3 k 6 + v 4 k s + . . . , (11) 
a 2 

which involves only even powers of k. The effective-range approximation corresponds to 
retaining only the first two terms in (11). In expansion (11), the quantities a and r are, 
respectively, the scattering length and the effective range, while the quantities v n are related 
to the potential shape. It should be noted that the quantity r is a measure of the effective 
interaction region. In the case of nucleon-nucleon interaction, the effective range r is ap- 
proximately equal to the range of nuclear forces (R). If, however, one deals with the doublet 
neutron-deuteron (nd) interaction, the effective range is anomalously large (tq ~ 500 fm), so 
that the statement that it is approximately equal to the range of nuclear forces (R ~ 2 fm) 
does not hold here. This is because, for nd interaction, the function k cot S has a pole in the 
vicinity of the point k 2 = owing to the existence of a low-energy virtual triton state [5]. 
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In a number of studies (see, for example, [6-8]), the authors approximated the S matrix 
by rational functions, this leading to so-called Bargmann potentials [4, 9]. The use of a 
Bargmann S matrix makes it possible to find explicit solutions to the direct and inverse 
scattering problems. Of particular importance are special cases where the Bargmann S matrix 
is determined by a small number of parameters such that physical observables are directly 
expressed in terms of these parameters. In [5], we showed that the use of the Bargmann S 
matrix corresponding to the presence of two states in the system being considered leads to 
the pole structure of the effective-range function (van Oers-Seagrave formula [10]). At low 
energies, this structure makes it possible to describe well doublet nd scattering on the basis 
of parameters that characterize the bound and the virtual state of the triton. 

In the present study, we use the Bargmann representation of the S matrix to construct 
the pole approximation of the effective-range function, bearing in mind that this pole ap- 
proximation is optimal for describing nucleon-nucleon scattering and that it involves a small 
number of parameters. The pole-approximation parameters have a clear physical meaning, 
and the values found for these parameters from an analysis of low-energy experimental data 
make it possible to obtain an excellent unified description of a bound state in the two-nucleon 
system (deuteron) and the triplet phase shift for neutron-proton scattering over a very wide 
energy range (0 — 1000 MeV). 

2. BARGMANN REPRESENTATION OF THE S MATRIX AND POLE 
APPROXIMATION OF THE EFFECTIVE-RANGE FUNCTION 

In order to investigate the interaction in the system of two particles, we will use the 
corresponding Bargmann S matrix [4, 9], which possesses simple analytic properties in the 
k plane. In this case, the Jost function F (k) can be chosen in the simplest way in the form 
of a rational function, 

/ j \ N k — ia n , . 

f (k = n - — r^. 12 

V ' n =l k + %\ n ' 

It has N simple zeros at the points k = ia n and N simple poles at the points k = —i\ n and 
exhibits the correct asymptotic behavior (4) at high energies. 

As a simple example of a Bargmann S matrix, one can consider the truncated S matrix 
corresponding to the well-known explicitly solvable Hulthen potential 

V(r) = -V {e r / R -iy 1 , (13) 



for which we have 



where 



g-n 2 n 
an = ^R> Xn= 2R> (14) 



g = -^V R 2 . (15) 

The Bargmann scattering matrix 

S(k) = U *±*»*±%, (16) 
n =i k — ia n k — i\ n 

which corresponds to the Jost function (12), has N simple "physical" poles at k = ia n in the 
k plane. The poles lying on the imaginary axis in the upper half-plane (a n > 0) correspond 
to bound states, while the poles in the lower half-plane correspond to resonances or virtual 
states. The resonances in the lower half-plane are grouped in pairs symmetric with respect to 
the imaginary axis, while the virtual states are on the imaginary axis (a n < 0). In addition to 
the physical poles corresponding to resonances and bound and virtual states, the Bargmann 
S matrix (16) has N so-called "redundant" poles at k = i\ n in the upper half-plane, which 
do not correspond to any physical states. The presence, along with physical poles, of the 
same number of redundant poles serves as some kind of a compensation factor, ensuring the 
correct asymptotic behavior [see Eq. (4)] of the Jost function at infinity. 

With allowance for Eqs. (2) and (3), the effective-range function (10) at real values of k 
can be expressed in terms of the real and the imaginary part of the Jost function as 

. it Re F (k) 

kcot8 = -k- — -7-4. (17) 
Im F (k) y ' 

From the Bargmann representation (16) of the S matrix, it immediately follows that the 
effective-range function kcotS can then be written in the form of a rational function of k 2 , 

kcot6 = TT^TwT ' ( 18 ) 

Q N -i {k 2 j 

where the degrees of the polynomials Pn ik 2 ) and Qn-i (k 2 ) of k 2 are N and N — 1, re- 
spectively, and where the coefficients in these polynomials are completely determined by the 
quantities a n and A n . 

The representation in (18) for the function k cot 5 can be considered as a Pade approxima- 
tion of this function [11] if the expansion of the function P/v ik 2 ) /Qn-i (k 2 ) = [N/ (N — 1)] 



in a Taylor-Maclaurin series coincides with the effective-range expansion (11) up to terms 
of order 2N — 1. In other words, the coefficients of 1, k 2 , ... , /c 2 ( 27V_1 ) i n the Taylor ex- 
pansion of the function [N/ (N — 1)] must coincide with the corresponding coefficients in 
the series in (11). The Pade approximation method was used by various authors to study 
nucleon-nucleon scattering (see, for example, [12-15]). It is of importance in our case that, 
from the Bargmann representation (16) of the S matrix, it automatically follows that the 
degree of the polynomial P in the numerator of the Pade approximant (18) must be greater 
by unity than the degree of the polynomial Q in its denominator, this being in accord with 
the condition of the theorem [12] on the solvability of the inverse scattering problem. In 
addition, one can see that, in the majority of cases, the condition L = M + 1 for the Pade 
approximant [L/M] of the function A; cot 5 is optimal for specific fits [13-15]. 

Let us now consider in detail the important particular case where the Bargmann S matrix 
corresponds to the presence of two physical states (N = 2) in the system. This is so, for 
example, in the doublet scattering of a neutron on a deuteron [5], in which case there are two 
triton states in the system, a bound and a virtual one. The scattering matrix corresponding 
to the presence of two states in the system can also be used to describe triplet neutron-proton 
scattering. Since the second state has not yet been observed experimentally in that case, the 
energy of this state must be considerably higher than the deuteron binding energy. If, in 
this case, the second state is a bound state of two nucleons, such a situation corresponds 
to the phenomenology of nodes that is described by a short-range deep attractive potential 
involving forbidden states [16-18]. 

If the system has two physical states, the rational Jost function and the Bargmann S 
matrix corresponding to it are given by 

F(*)= *-'°*-? , (19) 
k + tX k + tfi y ' 

k + ia k + i/3 k + i\ k + i\i 

k — ia k — i(3 k — i\ k — i\x 

The first and the second factors in the S'-matrix representation (20) correspond to either 

bound or virtual states of the system, while the third and the fourth factors correspond to 

redundant poles of the S matrix. The negative energies of the bound and virtual states of 

the system are 

h 2 a 2 , . 



h 2 B 2 

The expression describing the effective-range function k cot 5 and corresponding to the 
presence of two states in the system immediately follows from (17) and (19). The result is 

1 1 - c 2 k 2 + c 4 k 4 
kcot6 = - a 1 + Dk 2 > ( 23 ) 

where, for the sake of convenience, we have explicitly isolated the scattering length a. The 
parameters appearing in the pole representation (23) of the function k cot 5 and the param- 
eters a, (3, A, and \i of the Bargmann S matrix (20) are related by the equations 

1111 

a = - + - + - + -, 24 
a p X \i 

111111 

ap aX a>n pX p[i A/i 

c 4 = , (26) 

D = _l_a ± £ + X + v^ 
a apXfi 

The parameter D appearing in (23) determines the pole of the effective-range function k 2 : 

k l = -^- ( 28 ) 



Expression (23), which was obtained for the function A;cot<5 by the above method from 
the Bargmann representation of the S matrix, is a one-pole representation and involves four 
independent parameters. It can easily be shown that, apart from the form of presentation, the 
pole expression (23) coincides with the well-known empirical formulas that were given by van 
Oears and Seagrave [10] and by Cini, Fubini, and Stanghellini [19] and which are applied to 
describe, respectively, neutron-deuteron and nucleon-nucleon scattering. It is worth noting 
that various forms of the pole structure of the effective-range function k cot 5 have been 
repeatedly discussed for a rather long time and have been successfully used to describe 
neutron-deuteron [5, 10, 20-30] and nucleon-nucleon interactions [19, 31-35]. In the majority 
of cases, however, the formulas for k cot 5 were obtained empirically. But in our case, the 
pole formula for the function k cot 5 immediately follows from the Bargmann S matrix, 
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which possesses simple properties and which takes into account basic physical properties of 
the interaction in the system being considered. 

In [5], we showed that, in the case of neutron-deuteron interaction, the presence of the 
pole in the function k cot 8 is a direct consequence of the existence of a virtual triton state at 
a low energy. In the present study, the effective-range function having a pole structure will 
be used to describe neutron-proton interaction in the triplet (t) spin state 3 S±. In this case, 
the neutron-proton system has one bound state (deuteron), the scattering of a neutron on a 
proton at low energies (up to energies of about 10 MeV) being well described in the effective- 
range approximation. In this connection, it is convenient to recast the pole-approximation 
formula (23) into the form 

t crtf= __ + _ r ^ + __, (29) 

where a t and r t are, respectively, the scattering length and the effective range, while the 
parameter V2 determines the dimensionless shape parameter P t , which is widely used in the 
literature, according to the relation 

Pt = v 2 jr\. (30) 

We also note that, in the literature, definition of the dimensionless shape parameter P t 
frequently differs from the one in Eq. (30) by sign. But our definition here is more convenient 
for future discussion. 

The first two terms of the representation in (29) correspond to the effective-range ap- 
proximation, while the last, pole, term describes the deviation from this approximation. The 
presence of this pole term makes it possible to improve, with the aid of only two additional 
parameters, the description of the phase shift significantly and to extend the range of ap- 
plicability of this description greatly. Thus, the pole approximation (29), which was derived 
from the Bargmann representation of the S matrix, is a direct generalization of the effective- 
range approximation to the case where there are two physical states in the system. We note 
that the form (29) of the pole formula is more convenient for describing nucleon-nucleon 
scattering than that in (23), since, in (29), the low-energy scattering parameters a t , r t , and 
t>2 are isolated explicitly. 

The parameters a t and D appearing in (29) are related to the parameters of the Bargmann 
S matrix (20) by Eqs. (24) and (27), while the effective range r t and the shape parameter 
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i>2 are given by 

r t = -(D + c 2 ) , (31) 
a t 

»°=-0< + s)- (32) 

where c 2 and C4 are related to the ^-matrix parameters by Eqs. (25) and (26). 

Let us now consider the specific case where the S matrix in (20) for the case of two states 
reduces to the S matrix for one state, 

<*> 

In this case, the second state goes to infinity: (5 — > 00 and fj, — > 00. It can easily be seen 
that the coefficients D and v 2 then vanish in formula (29), which reduces, after that, to the 
effective-range-approximation formula 

k cot 5 = -— + -r t k 2 , (34) 
a t 2 

where the scattering length a t and the effective range r t are given by 

11 

<h = ~ + t , 35 
a A 

a + A 

The dimensionless asymptotic normalization factor Cd characterizing the bound state of 
the two-nucleon system (deuteron) can be expressed in terms of the residue of the S matrix 
at the pole k = ia as 

ZOi k=ta 

In the case of the S matrix for one state, Cd is given by 

03=^. ,3S, 

Upon expressing A in terms of a and Cj, formulas (35) and (36) for the scattering length a t 
and the effective range r t , respectively, can be recast into the form [36] 

2 C 2 
a 1 + Ci 



n = i(i-i). (40) 



Thus, formulas (39) and (40) give explicit expressions for the low-energy scattering param- 
eters a t and r t in terms of parameters that characterize the bound state of the two-nucleon 
system (deuteron) in the case where the interaction in the system is described by the S 
matrix for one state. The inverse statement is also valid — namely, the parameters charac- 
terizing the deuteron can be expressed, in this case, in terms of the scattering parameters a t 
and r t as 



1 

a = — 

n 



l 

at 

1 



2r t x 



(41) 



c* = m • ( 42 ) 

(l-2r t /a t ) V2 

The case where, at finite values of/? and /i, the coefficient D vanishes, while the parameter 
v 2 differs from zero is yet another important specific case of formula (29). This is so if the 
system has a virtual state whose wave number (3 is given by 

(3 = - ( a + A + fi) . (43) 

It can be seen from (27) that, in this case, the parameter D is equal to zero, while expression 
(29) for the effective-range function reduces to the expression corresponding to the shape- 
parameter approximation, 

k cot 5 = -— + -r t k 2 + v 2 k 4 , (44) 
a t 2 

where 

r t = 2- , (45) 
a t 

v 2 = Ptr\ = -- , (46) 
a t 

with a t , c 2 , and c 4 being given by (24)-(26). The effective-range approximation (34) and the 
shape-parameter approximation (44) are among the most popular and important methods 
for parametrizing data on nucleon-nucleon scattering at low energies [37, 38]. 

In the triplet spin state, the neutron-proton scattering length a t is positive. It follows 
that, in the approximation specified by Eq. (44), the sign of the shape parameter P t is 
determined by the sign of the parameter C4 in accordance with (46). Since the redundant 
poles k — iX and k = i\i of the S matrix lie in the upper half-plane and since a > and 
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(5 < 0, it follows from (26) that c 4 is negative. Therefore, the dimensionless shape parameter 
P t is positive in this approximation: 

Pt > . (47) 

Thus, we have shown that the shape-parameter approximation (44) directly follows from 
the Bargmann representation of the S matrix in the case where the system has two physical 
states of which at least one is virtual. If the system has two virtual states (a < 0, f3 < 0), 
the shape-parameter approximation (44) follows from the Bargmann representation of the 
S matrix under the condition 

(a + (3) = - (A + n) (48) 
and corresponds to the description of the neutron-proton interaction in the singlet spin state. 

3. CHOICE OF PARAMETERS AND DESCRIPTION OF SCATTERING IN THE 

TWO-NUCLEON SYSTEM 

An investigation of the phase shift S as function of energy plays a fundamental role in an 
analysis of data on nucleon-nucleon scattering. A traditional way to study this dependence 
at low energies consists in applying the effective-range approximation (34) and the shape- 
parameter approximation (44) [37, 38]. In order to determine the parameters a t , r t , and P t , 
one employs experimental data on scattering and also the experimental value of the deuteron 
binding energy e&. It was shown in [38] that, in this case, the determination of the scattering 
parameters a t , r t , and P t involved ambiguities, and this is at odds with the meaning of 
the effective-range expansion (11). The ambiguity in the determination of the low-energy 
scattering parameters was due primarily to an insufficient accuracy of experimental data 
at low energies. At the present time, the accuracy of experimental data is such that the 
ambiguity in the determination of the scattering length a t , the effective range r t , and the 
shape parameter P t is quite removable. 

By using the least squares method to construct a polynomial fit to the triplet function 
k cot 5 at low energies (T lab < 10 MeV) and relying on the latest experimental data of the 
VPI/GWU group of Arndt et al. on neutron-proton phase shifts [39], we obtain the following 
values of the triplet low-energy scattering parameters a t , r t , and v 2 : 

(H = 5.4030 fm, (49) 
11 



r t = 1.7494 fm, (50) 

v 2 = 0.163 fm 3 . (51) 

The resulting dimensionless shape parameter, P t — v 2 /rf — 0.0304, proved to be positive. 
This is in agreement with the estimate in (47), which was obtained on the basis of the 
Bargmann representation of the S matrix. 

In order to describe triplet scattering and the bound state in the two-nucleon system, 
we will use the pole approximation of the effective-range function (29). In doing this, we set 
the low-energy parameters a t , r t , and v 2 to the values in (49)-(51) and choose the value of 
the pole parameter D to be 

D = —0.225526 fm 2 . (52) 
With allowance for (28), this corresponds to the laboratory energy 

T = 368.026 MeV, (53) 

at which the experimental value of the triplet phase shift [39] changes sign. Thus, it can be 
seen, that, if the pole parameter D is fixed by using the experimental point 5 (T ) = 0, the 
quantity D can be considered as an experimental parameter directly determined to a rather 
high degree of accuracy. 

If one uses the parameter values in (49)-(52), then, as can be seen from the results of 
phase-shift calculations in Table 1 (the pole approximation PI), the pole approximation (29) 
describes the experimental triplet phase shift [39] up to laboratory energies of about 500 MeV 
with an absolute error not exceeding 0.5°. For the sake of comparison, the results obtained 
by calculating the phase shift within the effective-range approximation (ER) and within the 
shape-parameter approximation (SP) with the low-energy parameters set to the values in 
(49)-(51) are also presented in Table 1. As might have been expected, the effective-range 
approximation (34) describes well the phase shift only in the region of very low energies 
(Tab % 10 MeV). The shape-parameter approximation (44) describes well the experimental 
phase shift up to the energies of about 50 MeV (to within about 0.5°). At the same time, 
the pole approximation (29) provides an excellent description of the phase shift within a 
broad energy range, considerably improving the description in the low-energy region — for 
example, the experimental phase shift in the energy range T lab ^ 40 MeV is described by 
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formula (29) with an error not exceeding 0.01°, i.e. with better accuracy than accuracy of the 
experimental data (which is typically > 0.01°). In the above energy range, one can therefore 
consider the parametrization of the phase shift by the pole formula (29) as an alternative to 
the data of a partial- wave analysis (PWA). As a matter of fact, the approximation specified 
by Eq. (29) with the parameters given in (49)-(52) provides an excellent description of the 
phase shift in much wider energy interval than that which was used to find the parameters 
of this approximation. As can be seen from Fig. 1, the accuracy of the description of the 
phase shift undergoes virtually no deterioration up to an energy of 900 MeV, being 0.6° for 
laboratory energies of T lab ^ 900 MeV. Along with the phase shift corresponding to the pole 
approximation (29), the phase shift calculated within the effective-range approximation (34) 
and that calculated within the shape-parameter approximation (44) are also displayed in 
Fig. 1 for the sake of comparison. 

By slightly varying the shape parameter i> 2 and leaving the parameters a t , r t , and D 
unchanged, one can improve further the quality of description of the phase shift in the 
energy range T iab < 400 MeV. The minimum absolute error of description of the phase shift 
in this energy range is attained at 



and is about 0.1°, as can be seen in Table 1 (pole approximation P2). Thus, the pole ap- 
proximation (29) makes it possible to describe, by using a small number of parameters, the 
phase shift over a wide energy range to a precision close to that in determining experimental 
data. 

It should be noted that the value of the pole parameter D is very well determined by 
the experimental phase shift even at low energies and is close to the above value in (52). 
This can be shown most clearly by analyzing the dependence of the absolute error of the 
description of the phase shift, 



on the parameter D in a given range ^ T ^ T max . The dependence of the accuracy A on 
the parameter D for the energy interval ^ T ^ 20 MeV is displayed in Fig. 2. The phase 
shift in this case is described by the pole approximation (29) at fixed values of the low-energy 



v 2 = 0.168 fm 3 



(54) 



A =„ l<W ( T ) - <W CO | , 




(55) 
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parameters a t , r t , and v 2 (49)-(51), and at varying value of the pole parameter D. As can 
be seen from Fig. 2, the minimal error is achieved at negative values of D that are close to 
the value in (52), the dependence A (D) having quite a sharp character in the vicinity of 
the minimum, so that any significant deviation of the parameter D from the point of the 
minimum leads to a considerable deterioration of the quality of phase-shift description. We 
note that the value of D = corresponds to the shape-parameter approximation and that, 
for D — > oo, we obtain the effective-range approximation. The calculations show that, at 
the value of D = — 0.2147 fm 2 , which corresponds to the minimal deviation A, the absolute 
error of the description of the phase shift on the basis of (29) does not exceed 0.006° in the 
energy range being considered. At the same time, A takes values of 1.16° and 0.064° for, 
respectively, the effective-range and the shape-parameter approximation. Thus, we see that, 
even at low energies, the quality of the description of the phase shift by the pole formula is 
an order of magnitude higher than the quality of the description within the shape-parameter 
approximation, the value of the parameter D being negative and close to that in (52). We 
emphasize that the energy of T = 386.6 MeV, at which the calculated phase shift changes 
sign and which is in good agreement with the experimental value in (53), corresponds to the 
above value of D = — 0.2147 fm 2 . Calculated values of the pole parameters D and T , and of 
the accuracy A are given in Table 2 for different energy intervals. As can be seen from Table 
2, the pole parameter D is a quantity that admits an unambiguous determination yielding 
results that are rather weakly dependent on the interval of fitting. In fact, any fit to the 
phase shift leads, even at low energies, to negative values of the parameter D that are close 
to the value in (52). Thus, we see that, in just the same way as the parameters a t , r t , and 
i>2, D is a low-energy parameter that can be determined reliably and to a high precision. 

The pole formula (29) for the function k cot 5 coincides in form with the well-known 
empirical Cini-Fubini-Stanghellini (CFS) formula [19], which was used to describe nucleon- 
nucleon scattering [19, 31-35] and which is generally written as 



1 1 nh 

kcot6 = + r tk ?-JP. (56) 

a t 2 1 + qk 2 



A distinctive feature of the Cini-Fubini-Stanghellini approach is that the parameters p = —v 2 
and q = D in formula (56) are not independent — they are determined on the basis of the 
one-pion-exchange (OPE) theory of nucleon-nucleon interaction. Thus, the quantities p and 
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q in (56) are rather complicated functions of the scattering length a t , the effective range 
r t , the pion mass m w , and the pion-nucleon coupling constant G 2 . According to one-pion- 
exchange theory, the expressions for the parameters q and p can be written in the form [32, 
33] 

2 2-f 2 M (Z/V2-^/a t -r t /2K) 
9 K 1-PM(1/2V2-Xja t ) ' [b7) 

p = (A 2 - q/2) (2V2K -r t - 4A 2 /a t ) , (58) 

where M = /m n is the ratio of the nucleon mass to the pion mass, X w = h/m n c is the pion 
Compton wave length, and f 2 = (m n /2m N ) 2 G 2 . For the low-energy scattering parameters 
a t and r t , we will use, in the following, the values obtained above and quoted in (49) and 
(50). On this basis, one can readily calculate the shape parameter i> 2 and the pole parameter 
D within the Cini-Fubini-Stanghellini approach. The results are 

y 2 CFS = -0.121 fm 3 , (59) 

D CFS = 3.777 fm 2 . (60) 

It can be seen that the parameters V2 and D calculated within the Cini-Fubini- 
Stanghellini approach differ considerably from the "experimental" values that are quoted 
in (51) and (52) and which are determined quite reliably from present-day data on the 
triplet phase shift [39]. It should be noted that the distinction is not only quantitative but 
also qualitative since the parameters v 2 and D have opposite signs within the Cini-Fubini- 
Stanghellini approach. The phase shift is described poorly with the Cini-Fubini-Stanghellini 
parameters. The explanation for so sharp a discrepancy with experimental data is likely to be 
the following. It has been firmly established that one-pion exchange is not the only mecha- 
nism and even is not the main mechanism of nucleon-nucleon interaction — the contribution 
of other mechanisms to nucleon-nucleon interaction is much more significant. Therefore, it 
comes as no surprise that the oversimplified one-pion-exchange scheme predicts erroneous 
values for the parameters V2 and D. For want of a theory that would make it possible to 
calculate the low-energy parameters of nucleon-nucleon interaction on the basis of a micro- 
scopic approach (QCD), one has to treat them as adjustable parameters that are determined 
directly from experimental data. 
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Along with the aforesaid, we note that the recent calculation of the shape parame- 
ter v 2 within effective field theory (EFT) in [40, 41] also yielded an incorrect sign of v 2 : 
V VFT _ —0.95 fm 3 , although calculations within this theory generally led to good agreement 
with experimental data for many computed features of the nucleon-nucleon and neutron- 
deuteron systems (see references quoted in [40, 41]). These discrepancies indicate that the 
shape parameter i> 2 and the higher order parameters v n are rather subtle and sensitive char- 
acteristics of the nucleon-nucleon interaction. In [40, 41], it was also indicated that the 
contribution to the nucleon-nucleon interaction from more "short-range" mechanisms than 
the one-pion-exchange mechanism is of importance. 



4. LOW-ENERGY PARAMETERS OF THE EFFECTIVE-RANGE EXPANSION 



Much attention has permanently been given to studying the low-energy parameters of 
the effective-range expansion (11) [12-15, 32-35, 37, 38, 40-50]. It should be noted that, 
while the scattering length a and the effective range r can be determined directly from 
experimental data to a fairly high degree of precision [42-45], the higher order parameters 
v n (n = 2, 3, 4, . . .) are less convenient for an experimental determination, their theoretical 
calculation becoming more involved as the parameter order increases. At the same time, the 
shape parameters v n , along with the scattering length a and the effective range r , are of 
particular importance for constructing and comparing various realistic models of nucleon- 
nucleon interaction and for describing nucleon-nucleon scattering. Moreover, the shape pa- 
rameters are of importance for investigating the physical properties of the bound state in 
the two-nucleon system (deuteron) — in particular, their role in the useful expansion for 
the root-mean-square radius of the deuteron was demonstrated in [51-55]. The correlations 
between the properties of the deuteron and the parameters of nucleon-nucleon scattering 
that are determined by the expansion in (11) are also studied [51-56]. It should be not- 
ed that, in the past years, the parameters v n were discussed and calculated in considering 
nucleon-nucleon interaction within effective field theory [40, 41]. 

In connection with the aforesaid, the determination of the parameters of the effective- 
range expansion (11) is of great importance. An explicit expansion of the function kcot5 in 
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a power series in k 2 follows directly from the pole representation (29): 

oo 

k cot 5 = -- + -r t k 2 + V v 2 {-Df- 2 k 2n . (61) 

at 2 n=2 

Therefore, all of the shape parameters v n can be also obtained explicitly for this case. We 
have 

v n = (-l) n v 2 D n ' 2 , n^3. (62) 

Thus, we see that, in the pole approximation, all parameters v n (n = 3, 4, 5, . . .) of the 
effective-range expansion are determined explicitly in terms of the shape parameter v 2 and 
the pole parameter D according to the simple formula (62) and can easily be calculated to 
any order n. 

It should be noted that the dimensionless shape parameters p n are also used along with 
the dimensional shape parameters v n [see Eq. (30)]. We define the dimensionless shape 
parameters p n according to the relation 

v n = p n r 2n - 1 . (63) 

It can easily be seen that the parameters p n are coefficients of the power expansion of the 
function kr t cot 5 t with respect to the dimensionless parameter x = kr t . The dimensionless 
shape parameters p n are convenient for qualitative analysis of the dependence of the phase 
shift on energy. In our case, the parameter D is negative by virtue of (52). Since we also 
have in our case p 2 = P t = 0.0304 > and \D/r 2 \ < 1, one can easily see that all of the 
parameters p n , defined according to Eqs. (62) and (63), are positive and decrease in absolute 
value; that is, 

p 2 > p 3 > ■ ■ ■ > Pn > Pn+1 > ■■■ > . (64) 

In Table 3, the low-energy scattering parameters a t , r t , and v 2 are presented along with 
the parameters v 3 , t> 4 and v 5 calculated by formula (63) at the experimental value of D = 
— 0.225526 fm 2 and the shape-parameter values of v 2 = 0.163 fm 3 (approximation PI) and 
v 2 = 0.168fm 3 (approximation P2). For the sake of comparison, the parameters of the 
expansion in (11) that were found in [50] by using the partial- wave-analysis data on nucleon- 
nucleon scattering that were obtained by the Nijmegen group [57] (Nijm version) are also 
given in Table 3. We note that, at the present time, the partial-wave-analysis data obtained 
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by the VPI/GWU group of Arndt et al. [39] and the Nijmegen group [57] are the most 
accurate and frequently used data on the phase shifts for nucleon-nucleon scattering. It can 
be seen from Table 3 that the values found here for the shape parameters v n (n = 3, 4, 5) 
on the basis of the pole approximation of the effective-range function are rather stable 
with respect to the variations in the shape parameter v 2 . It can easily be shown that the 
parameters v n are also rather weakly sensitive to the variation in the pole parameter D. 
Thus, formula (62) ensures a stable determination of the parameters v n . 

The following important comment concerning the applicability and accuracy of the simple 
one-pole approximation (29) for determining the parameters v n is in order. Since the smooth 
interpolation curve specified by (29) provides an excellent description of experimental data of 
Arndt et al. on the phase shift [39] over a wide energy range (in particular, to a precision not 
poorer than 0.01° for energies up to T lab = 40MeV), this curve determines the parameters 
v n with the degree of reliability and stability as high as that to which they can in principle 
be determined by present-day experimental data on the phase shift, the error in these data 
being > 0.01°. This means that, if a different interpolation curve also described well the 
experimental phase shift and, for some parameters v n , gave values considerably differing 
from those obtained here, one could conclude that these parameters v n are not determined 
by experimental data. As a matter of fact, the calculations show that only the parameters 
to V4 inclusively are determined more or less reliably. 

Table 3 shows that the results for the low-energy scattering parameters, obtained here 
for the partial-wave-analysis data presented by Arndt et al. [39] are radically different from 
the results of the calculations performed in [50] for the partial-wave-analysis data of the 
Nijmegen group [57], the distinction being not only quantitative but also qualitative — that 
is, all shape parameters calculated here on the basis of the phase shifts of Arndt et al. are 
positive and decrease, while the analogous parameters calculated for the phase shifts of the 
Nijmegen group increase in absolute value and include a negative parameter (t> 4 ). To be 
more exact, the latter observation concerns the dimensionless shape parameters p n , which 
can be easily calculated according to Eq. (63) and which are given in Table 4. 

A detailed investigation of the parameters of the effective-range expansion (11) will 
be given in our future publications. Here, we restrict ourselves to mentioning briefly the 
following: in order to reveal possible reasons behind such a discrepancy between the results, 
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we have studied the stability of the calculation of the shape parameters v n to variations in 
the triplet scattering length a t . We have shown that their sensitivity to variations in this 
quantity is extremely high. In particular, a change in the scattering length a t as small as a 
few tenths of a percent can lead to a severalfold change in the shape parameter v 2 . This is in 
accord with the comment in the preceding section that the shape parameter v 2 , as well as the 
higher order parameters v n , is a rather subtle and sensitive characteristic of nucleon-nucleon 
interaction. Thus, so sharp a distinction between the shape parameters v n for the GWU and 
the Nijmegen phase shifts is due to a significant difference in the scattering length: the value 
of af wlJ = 5.403 fm, which was obtained in the present study, differs from af ljm = 5.420 fm 
more than by 0.3%. This difference in the values of the scattering length leads to a decrease 
in the shape parameter v 2 by a factor of 4 for the phase shifts of the Nijmegen group in 
relation to the phase shifts of the GWU group — that is, from the value of w^ wu = 0.163 fm 3 
to the value of w^ ljm = 0.040 fm 3 . The discrepancy between the corresponding higher order 
parameters v n is still larger. At the same time, Table 3 shows that the values of the effective 
range r t for the GWU and the Nijmegen phase shifts are close to each other. It should also 
be noted that our results on the shape parameter v 2 are close to the results of some earlier 
studies. In particular, the value of the shape parameter v 2 = 0.137fm 3 corresponding to the 
dimensionless shape parameter P t = 0.027 for the Reid soft-core potential RSC [58] agrees 
well with the value of v 2 = 0.163fm 3 , which was obtained in the present study. 

We note that the Nijmegen value of the triplet scattering length af ljm = 5.420 fm is 
rather close to the presently recommended experimental value [59] 

af pt = 5.424 fm, (65) 

while the value calculated here on the basis of the GWU phase shifts [39], af WXJ = 5.403 fm, 
is close to some of the experimental values obtained previously for the triplet scattering 
length [32, 42, 43]. It should be emphasized that, for the experimental values of the triplet 
scattering length a t , various authors [32-35, 38, 39, 42-46, 58-61] present values changing 
within a broad range — from 5.37 [38, 60] to 5.479 fm [61]. We also note that the present-day 
experimental value of the triplet scattering length in (65) leads to exaggerated (in relation 
to experimental data) values of the root-mean-square radius of the deuteron [52, 53, 
56] and the asymptotic normalization constant As for it [56]. From all of the aforesaid, 
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we can therefore draw the following important conclusion: in order to remove the existing 
discrepancies between the current experimental value of a t in (65), on one hand, and the 
present-day values of the quantities and As for the deuteron and the partial-wave-analysis 
data of the GWU group, which lead to the value of af wu = 5.403 fm, on the other hand, it 
is of paramount importance to refine the experimental value of the triplet scattering length 
a t . The value of the triplet scattering length a t is also of particular importance since it is 
often used as one of the input values in fitting the parameters of nucleon-nucleon potentials. 
We note that the value of af WXJ = 5.403 fm, which corresponds to the partial-wave analysis 
performed by Arndt et al., is in perfect agreement with the experimental values of the 
quantities and As for the deuteron [56]. 

In a wide energy range [approximately to the energy corresponding to the pole of the func- 
tion kcotS (53)], the experimental triplet phase shift [39] is well described by the effective- 
range expansion (11) with the low-energy scattering parameters a t , r t , v 2 , v 3 , . . . set to the 
values determined here. This is a direct corollary of the fact that, for the phase shift borrowed 
from [39], all of the dimensionless shape parameters p n appear to be decreasing in absolute 
value. For example, a fifth-degree polynomial in k 2 describes the experimental phase shift in 
the energy range T\ ah ^ 200 MeV to within about 1° and in the energy range Ti ab ^ 50 MeV 
to within about 0.02°. The description of the function kcotS by various-degree polynomials 
corresponding to various polynomial approximations of the effective-range expansion (11) 
are shown in Fig. 3. We see that a successive increase in the degree of a polynomial leads to 
an improved description of the function k cot 5 (and, accordingly, of the phase shift) and to 
the extension of the interval where this description is valid. Thus, the use of the effective- 
range expansion with a small number of terms provides a good description of the phase shift 
over a wide energy range (0 — 250 MeV) if the low-energy parameters are set to the values 
found here. This indicates that, in contrast to statements advocated in some articles (see, for 
example, [50]), the effective-range expansion is highly accurate and very useful in describing 
nucleon-nucleon scattering. 

The potential of the effective-range expansion (11) is directly related to the radius of 
its convergence. The function k cot 5 considered in the complex plane of k 2 is an analytic 
function of k 2 within some region near the origin of coordinates; therefore, it can be expanded 
in a Taylor-Maclaurin series (11) in powers of k 2 in the vicinity of the point k 2 = 0. Thus, a 
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particular importance of the function k cot 5 is associated with its analyticity near the point 
k 2 = 0. Noyes and Wong [62] showed that the one-pion-exchange model for nucleon-nucleon 
interaction leads to the appearance of a cut in the scattering amplitude on the negative 
axis of energy for k 2 ^ —m 2 c 2 /4:h 2 ; therefore, the radius of convergence of the series in 
(11) is extremely small in the one-pion-exchange model, T = m 2 /2m,N = 9.7 MeV. This 
contradicts the fact that the experimental phase shift [39] is well described by the pole 
formula (29) with the parameters set to the values in (49)-(52), since, for the radius of 
convergence of the series in (11), formula (29) gives the value of T = 368.026 MeV, which 
we chose on the basis of the "experimental" condition S (T ) = (the convergence circle is 
determined by the condition \k 2 \ < 1/ \ D\). In connection with this contradiction, we recall 
the comments at the end of the preceding section that concern the limited applicability of the 
one-pion-exchange mechanism. As was indicated by Noyes himself [33], there is no a priori 
method for estimating the accuracy or the region of applicability of the phenomenological 
expansion (11), so that this question must be solved on the basis of experimental data. 

Numerical values of the dimensionless shape parameters p n provide direct information 
about the radius of convergence of the series in (11). In the case under consideration, the 
increase in the parameters p n in absolute value indicates, as is the case for the triplet np 
phase shift of the Nijmegen group, that the effective-range expansion (11) has a small radius 
of convergence, so that it is not very useful [50]. On the contrary, the decrease in the shape 
parameters p n for the triplet np phase shift of the GWU group suggests that the radius 
of convergence of the series in (11) is large and that this expansion is very useful in this 
case. This is corroborated by a high quality of the description of the phase shift in the case 
where the function k cot 5 is approximated by polynomials of relatively low degrees. Thus, 
the present-day data from the partial-wave analyses performed by the two main groups 
in [39] and [57] do not agree with each other both in what is concerned with the low- 
energy scattering parameters calculated on their basis and in what is concerned with the 
applicability of the effective-range expansion to them. At the same time, it is important to 
note that the numerical values obtained for the phase shifts by the GWU and the Nijmegen 
group are rather close to each other. This indicates once again that the low-energy scattering 
parameters are subtle and sensitive characteristics of nucleon-nucleon scattering. 
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5. DESCRIPTION OF THE PROPERTIES OF THE DEUTERON 



Following the same line of reasoning as in the case of constructing the effective-range 
expansion (11), one can write an expansion of the function A; cot 5 in a power series at the 
point k 2 = —a 2 — that is, at the energy equal to the deuteron binding energy e d = h 2 a 2 /rriN ■ 
This expansion has the form 

k cot 5 = -a + \-p d (k 2 + a 2 ) + w 2 (k 2 + a 2 ) 2 + ... , (66) 

where Pd = P (— ~~ e d) is the deuteron effective range corresponding to S'-wave interaction. 
The definition and the properties of the effective deuteron range pd and of the function 
p(Ei,E 2 ) are discussed in detail elsewhere [38]. Using Eqs. (29) and (66), we can easily 
establish that the quantities a and p d for the deuteron are related to the parameters of the 
pole representation of the effective-range function as 

1 , 1 2 U 2« 4 ,_-v 

Of 2 1 — Da 2 

2v 2 a 2 

P d = P m ~ n . ( 68 ) 

(1 — Z9a 2 ) 

where p m = p (0, — e d ) is the so-called mixed effective range given by 

Pm = -(l- —) ■ (69) 



a \ aa t/ 

Using formulas (67) and (69) , we can recast the expression for the mixed effective range p m 
into the form 

p m = r t - T -^. (70) 

From formulas (70) and (68), one can easily obtain expansions of the standard scattering 
effective range r t and the deuteron effective range pd in power series in a 2 . We have 

oo 

r, = p m -2^(-l)%„ +1 a 2 ", (71) 

71=1 

OO 

Pd = p m + 2J2 (-l) n nv n+1 a 2n , (72) 

n=l 

where the parameters v n are given by (62). The dimensionless asymptotic normalization 
factor Cd for the deuteron is expressed in terms of the deuteron effective range pd as 

Cj = (1 - ap d )- 1 . (73) 
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The constant Cd and the asymptotic normalization factor A$, which is widely used in the 
literature [47, 63], are related by the equation 

A% = 2aCl . (74) 

For various approximations of the effective-range function, we have calculated the fol- 
lowing parameters characterizing the deuteron: the binding energy Ed, the effective range pd, 
and the asymptotic normalization factors Cd and As- In Table 5, the results of these cal- 
culations are given along with their experimental counterparts from [64, 65]. We note that 
the deuteron binding energy was calculated by using a relativistic formula, which is more 
accurate than (21). It can be seen from Table 5 that the features of the deuteron that were 
calculated in the pole approximation (29) with the parameter values from (49)-(52) are in 
good agreement with their experimental values. The results obtained in the shape-parameter 
approximation (SP) differ insignificantly from their counterparts in the pole approximation, 
while the results of the calculations that take into account the cubic term in energy in the 
effective-range expansion are nearly coincident with those in the pole approximation. It can 
be seen from Table 5 that the convergence of the calculated features of the deuteron versus 
the number of terms that are taken into account in the effective-range expansion is very fast 
— the shape-parameter approximation yields a highly precise result, the inclusion of higher 
order terms in energy introduces virtually no changes in this result. This is a consequence of 
the fact that the dimensionless shape parameters p n are decreasing quantities in this case. 
It should be emphasized that the features of the deuteron were calculated on the basis of 
the scattering-parameter values in (49)-(52), which correspond to the phase shifts of Arndt 
et al., and this means that we have very good agreement between the experimental data of 
the GWU group on scattering [39] and the experimental data on the bound state (deuteron) 
from [64, 65]. 

The values found for the features of the deuteron in the effective-range approximation 
(ER) are somewhat exaggerated in relation to their experimental counterparts and the re- 
sults obtained in the pole approximation. This is due primarily to an insufficiently accurate 
determination of the deuteron binding energy in the effective-range approximation with the 
low-energy parameters a t = 5.403 fm and r t = 1.7494 fm. If one uses the experimental values 
for the deuteron binding energy and for the triplet scattering length (ed = 2.224589 MeV 
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and a t = 5.403 fm, respectively), the effective-range approximation yields the following val- 
ues for the deuteron effective range pd and for the asymptotic normalization factor Ag\ 
Pd = 1.7331 fm and As = 0.8795 fm -1 / 2 , these results being in good agreement with the 
corresponding experimental values and with the results of the calculations in the pole ap- 
proximation. This is in accord with the results obtained in [56], where it was established 
that the asymptotic normalization factor As depends only slightly on the model of nucleon- 
nucleon interaction; at the experimental value of the deuteron binding energy Ed, 99.7% of 
it is determined by the triplet scattering length a t . 



The parameters a, (3, A, and p of the Bargmann S matrix (20), which corresponds to the 
presence of two physical states in the system, are unambiguously related to the parameters 
of the pole approximation of the effective-range function (29) and, as can easily be seen, are 
the roots of the fourth-degree algebraic equation 



Solving Eq. (75) with the parameters set to the values in (49)-(52), we obtain the following 
values for the parameters of the Bargmann S matrix: 

a = 0.2315 fm" 1 , (3 = 1.2293 fm" 1 , A = 2.5603 + i3.5248 fm -1 , p = 2.5603 - i3.5248 fm -1 . 



In the case being considered, the two-nucleon system has two bound states. In accordance 
with the phenomenology of nodes that employs potentials involving forbidden states [16-18], 
the lowest, deeply lying, state characterized by the energy 
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(75) 



(76) 



£o = 63.7555 MeV, 



(77) 



is unobservable, while the excited state of binding energy 



e d = 2.2237 MeV 



(78) 



corresponds to the deuteron. 
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Thus, the use of the pole approximation in describing triplet neutron-proton scattering 
automatically leads to considering deep potentials involving forbidden states. 

In the shape-parameter approximation (44) (D = 0), we have the following values for 
the parameters of the Bargmann S matrix (20): 

a = 0.2315 fin -1 , 13 = -2.7787 far 1 , A = 1.2736 + i0.3784fm _1 , fi = 1.2736 - i0.3784 far 1 . 

(79) 

In this case, the second state is a virtual state at the energy 

e v = 353.4732 MeV, (80) 

while the ground state, whose binding energy is 

Ed = 2.2236 MeV (81) 

[in fact, it coincides with that in (78)], corresponds to the deuteron. As might have been 
expected, the parameters in (79) satisfy relation (43). 

7. CONCLUSION 

Our basic results and conclusions can be formulated as follows. Relying on the Bargmann 
representation of the S matrix, we have formulated the pole approximation for the effective- 
range function kcotS. At specific values of the S'-matrix parameters, the effective-range 
approximation and the shape-parameter approximation immediately follow from this approx- 
imation. The pole approximation of the function k cot 8 is optimal for describing nucleon- 
nucleon scattering and involves a few parameters. The parameters of this approximation 
have a clear physical meaning. They are related to the parameters of the Bargmann S ma- 
trix by simple equations. It has been shown that the pole approximation deduced from the 
Bargmann representation of the S matrix is a direct generalization of the effective-range 
approximation to the case where there are two physical states in the system. The presence 
of the pole term makes it possible to improve significantly, by using only two additional 
parameters, the description of the phase shift and to expand sizably the applicability range 
of this description. In the shape-parameter approximation corresponding to the Bargmann 
representation of the S matrix, one can obtain an important constraint on the dimensionless 
shape parameter, P > — that is, the shape parameter is positive. 
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By using a least squares polynomial fit to the function k cot S at low energies, we have 
obtained, on the basis of the analysis of the latest experimental data on phase shifts that 
was performed by the VPI/GWU group, the triplet low-energy parameters of neutron-proton 
scattering: a t = 5.4030 fm, r t = 1.7494 fm, and v 2 = 0.163 fm 3 . With these values of a t , r t , 
and v 2 and the pole parameter D, the pole approximation of the function A; cot 5 provides 
an excellent description of the triplet phase shift for neutron-proton scattering over a wide 
energy range (Ti ab < 1000 MeV), the description in the low-energy region also being improved 
considerably. 

For the experimental phase shifts of the VPI/GWU group, the values that we have 
obtained on the basis of the pole approximation for the triplet shape parameters of the 
effective-range expansion are positive and decrease with increasing n. On the basis of the 
effective-range expansion with the values found for the low-energy scattering parameters a t , 
r t, v 2 , vz, ■ ■ ■ , the phase shift is described well over a wide energy range extending approxi- 
mately to the energy at which the phase shift changes sign, this being a direct consequence 
of a decrease in the dimensionless shape parameters with increasing n. This circumstance 
is indicative of a high precision of the effective-range expansion and its high potential for 
describing experimental data on nucleon-nucleon scattering, in contrast to the statements 
of some authors (see, for example, [50]). 

The results obtained here for the shape parameters by using the data of the partial- wave 
analysis performed by the GWU group differ drastically from the results of the calculations 
in [50] for the data of the partial-wave analysis performed by the Nijmegen group, this 
distinction being not only quantitative but also qualitative — that is, all of the dimensionless 
shape parameters calculated here by using the phase shifts of the GWU group are positive and 
decrease, while the analogous parameters calculated for the Nijmegen phase shifts increase in 
absolute value and include a negative parameter (pa). In our opinion, so sharp a discrepancy 
between the shape parameters for the GWU and the Nijmegen phase shifts is due to quite a 
significant difference in the scattering length: the value af wu = 5.403 fm obtained here differs 
from af ^ m = 5.420 fm more than by 0.3%. This difference in the scattering length leads to the 
decrease in the shape parameter v 2 by a factor of 4 for the Nijmegen phase shifts in relation 
to the phase shifts obtained by the GWU group — that is, from the value u^ wu = 0.163 fm 3 
to the value of v 2 ^ m = 0.040 fm 3 — the discrepancy between the corresponding values of 
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the higher order parameters v n being still greater. This confirms that the shape parameter 
i>2, as well as the higher order parameters v n , is a rather subtle and sensitive characteristic 
of nucleon-nucleon interaction. The aforesaid leads to the important conclusion that an 
experimental refinement of the triplet scattering length a t is of paramount importance since 
it is necessary to remove the existing discrepancies between the current experimental value 
of af pt = 5.424 fm, on one hand, and the present-day values of and As for the deuteron 
[56] and the present-day data of the partial-wave analysis performed by the GWU group, 
which lead to the value of ap wu = 5.403 fm, on the other hand. 

For various approximations of the effective-range function, we have calculated the main 
features of the deuteron — the binding energy Ed, the effective range pd, and the asymptotic 
normalization factors Cd and A$. The results obtained for them in the pole approximation 
agree very well with their experimental counterparts. The results in the shape-parameter ap- 
proximation differ insignificantly from those in the pole approximation. We have found that 
the convergence of the calculated features of the deuteron versus the number of terms that 
are taken into account in the effective-range expansion is very fast — the shape-parameter 
approximation gives a nearly precise result, which undergoes virtually no changes upon tak- 
ing into account higher order terms in energy. Thus, we can state that, on the basis of the 
Bargmann representation of the S matrix, a good unified description has been obtained for 
the bound state of the two-nucleon system and the triplet phase shift for neutron-proton 
scattering up to energies of about 1000 MeV. 
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Table 1. Triplet phase shift for np scattering as a function of the laboratory energy T lab ac- 
cording to calculations within the effective-range (ER) approximation, the shape-parameter 
(SP) approximation, and the pole (PI and P2) approximations (for the shape-parameter 
values v 2 = 0.163 fm 3 and v 2 = 0.168 fm 3 , respectively). 



MeV 


Phase shift 5 t , deg 


ER 


SP 


PI 


P2 


Experiment [39] 


1 


147.84 


147.83 


147.83 


147.83 


147.83 


5 


118.34 


118.23 


118.23 


118.23 


118.23 


10 


102.93 


102.56 


102.55 


102.54 


102.55 


25 


81.87 


80.36 


80.26 


80.20 


80.26 


40 


71.19 


68.44 


68.11 


68.01 


68.11 


50 


66.23 


62.68 


62.14 


62.01 


62.11 


75 


57.51 


52.16 


50.90 


50.70 


50.77 


100 


51.64 


44.80 


42.60 


42.34 


42.34 


125 


47.30 


39.25 


35.95 


35.66 


35.58 


150 


43.93 


34.89 


30.37 


30.05 


29.94 


175 


41.19 


31.35 


25.53 


25.19 


25.09 


200 


38.92 


28.42 


21.23 


20.90 


20.84 


225 


36.99 


25.95 


17.36 


17.04 


17.03 


250 


35.32 


23.84 


13.81 


13.53 


13.57 


275 


33.86 


22.01 


10.53 


10.29 


10.37 


300 


32.57 


20.42 


7.47 


7.28 


7.38 


325 


31.41 


19.02 


4.60 


4.47 


4.56 


350 


30.38 


17.78 


1.88 


1.82 


1.87 


375 


29.44 


16.68 


-0.71 


-0.69 


-0.71 


400 


28.58 


15.69 


-3.18 


-3.07 


-3.20 


425 


27.79 


14.79 


-5.54 


-5.34 


-5.61 


450 


27.06 


13.99 


-7.82 


-7.51 


-7.94 


475 


26.39 


13.25 


-10.02 


-9.59 


-10.21 


500 


25.77 


12.58 


-12.14 


-11.60 


-12.41 



32 



Table 2. Calculated values of the pole parameters D and T , and of the accuracy of the 
phase-shift description A for different energy intervals ^ T ^ T max . 



T max , 


D, 

7 


T , 

yj ? 


A, 


MeV 


fm 2 


MeV 


dee 

o 


10 


-0.2007 


413.5 


0.0056 


20 


-0.2147 


386.6 


0.0057 


30 


-0.2172 


382.1 


0.0058 


40 


-0.2236 


371.2 


0.0091 


50 


-0.2303 


360.4 


0.014 


60 


-0.2355 


352.4 


0.020 


80 


-0.2416 


343.5 


0.029 


100 


-0.2432 


341.3 


0.032 


150 


-0.2412 


344.1 


0.061 


200 


-0.2356 


352.3 


0.16 


250 


-0.2315 


358.5 


0.26 


300 


-0.2290 


362.4 


0.33 


350 


-0.2279 


364.2 


0.36 


368 


-0.2256 


368.0 


0.44 
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Table 3. Triplet low-energy scattering parameters obtained by using the data of the partial- 
wave analysis performed by the GWU group (versions GWU PI and GWU P2 correspond to 
the pole approximation at the shape-parameter values of v-i = 0.163fm 3 and v 2 = 0.168 fm 3 , 
respectively) and the data of the partial-wave analysis performed by the Nijmegen group 
[50, 57] (Nijm version). 



Version 


a t , fm 


r u fm 


i>2, fm 3 


t>3, fm 5 


i>4, fm 7 


i>5, fm 9 


GWU PI 


5.4030 


1.7494 


0.163 


0.037 


0.0083 


0.0019 


GWU P2 


5.4030 


1.7494 


0.168 


0.038 


0.0085 


0.0019 


Nijm 


5.420 


1.753 


0.040 


0.672 


-3.96 


27.1 
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Table 4. Triplet dimensionless shape parameters p n obtained by using the data of the 
partial-wave analysis performed by the GWU group (version GWU corresponds to the pole 
approximation PI at the shape-parameter value of t>2 = 0.163 fm 3 ) and the data of the 
partial- wave analysis performed by the Nijmegen group [50, 57] (Nijm version). 



Version 


P2 


P3 


P4 


P5 


GWU 


0.0304 


0.0023 


0.00017 


0.000012 


Nijm 


0.00743 


0.0406 


-0.0778 


0.173 
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Table 5. Features of the deuteron according to calculations within various approximations 
of the effective-range expansion (ER, SP, v n ) and in the pole approximation PI. 



Version 


e r u MeV 


Oru fm 




A s , fm- 1 / 2 


ER 


2.2387 


1.7494 


1.2979 


0.8846 


SP 


2.2236 


1.7145 


1.2876 


0.8761 


^ o 


2.2237 


1.7151 


1.2878 


0.8763 


1>4 


2.2237 


1.7151 


1.2878 


0.8763 


V5 


2.2237 


1.7151 


1.2878 


0.8763 


PI 


2.2237 


1.7151 


1.2878 


0.8763 


Experiment 


2.224589 


1.7251 


1.2904 


0.8781 
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Fig. 1. Triplet phase shift for neutron-proton scattering as a function of the laboratory 
energy according to calculations within the effective-range (ER), the shape-parameter (SP), 
and the pole (PI) approximation. The points represent experimental data borrowed from 
[39]. 
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Fig. 2. Dependence of the accuracy of the phase-shift description A on the pole parameter 
D for the energy interval < T < 20 MeV . 
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Fig. 3. Triplet effective-range function kcotS versus the laboratory energy for various poly- 
nomial approximations of the effective-range expansion (11). The figures on the curves indi- 
cate the order of approximation. The upper curve corresponds to experimental data borrowed 
from [39]. 
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